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PREFACE 

This  handbook  presents  the  methods  used  in  testing  an  inertial  navi- 
gation system  at  the  Air  Force  Flight  Test  Center  (AFFTC) , Edwards  AFB, 
California.  The  work  was  done  under  the  authority  of  the  Improved 
Navigation  Systems  Testing  and  Analysis  Study  Plan. 

The  format  of  this  handbook  was  chosen  to  make  it  usable  to  project 
engineers  of  the  Systems  Engineering  Branch  at  the  AFFTC.  As  such, 
information  is  presented  to  give  a novice  in  the  field  of  inertial 
navigation  system  evaluation  sufficient  background  and  knowledge  to 
perform  an  accurate  evaluation  of  inertial  navigation  systems. 

The  authors  wish  to  acknowledge  the  following  individuals  who  were 
instrumental  in  the  preparation  of  this  handbook:  Mr  B.  Lyle  Schofield, 

Chief,  Flight  Test  Technology  Branch,  for  guidance  and  editorial  comments, 
and  Mr  William  Taylor,  Mathematician,  for  assisting  in  the  development 
of  the  computer  software. 
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INTRODUCTION 


The  purpose  of  this  document  is  to  provide  the  Flight 
Test  Engineer  with  a procedural  document  for  the  planning 
and  conduct  of  a flight  test  evaluation  of  any  inertial 
navigation  system  (INS) . Included  in  this  document  is  a 
description  of  the  types  of  inertial  navigators  that  will 
be  tested  at  the  Flight  Test  Center  in  the  foreseeable 
future  and  guidelines  on  how  to  plan  the  testing  from  an 
initial  estimate  of  the  required  flying  hours  to  the  writing 
of  the  final  report.  The  detailed  test  procedures  required 
to  collect  the  necessary  data  are  covered  including  sample 
flight  cards.  Data  collection  and  analysis  are  discussed 
using  actual  test  data  to  substantiate  the  techniques 
addressed.  Methods  of  presenting  the  data  in  various  re- 
ports are  discussed. 

The  appendixes  contain  information  on  the  two  computer 
programs  used  to  analyze  the  test  data,  NAVAN  and  CEPLOT. 

The  program  description  and  users  guide  is  contained  in 
appendixes  A and  B,  with  check  cases  in  appendix  C. 


ORY  OF  OPERATION 


Inertial  space  is  defined  as  that  space  where  Newton's 
laws  of  motion  apply.  Inertial  navigation  is  based  upon 
measurements  made  with  respect  to  inertial  space.  An 
inertial  system^determ ines  the  displacement  of  the  carrying 
vehicle  from  its  starting  point  by  measuring  the  accelerations 
of  the  vehicle  relative  to  the  earth  and  integrating  the 
accelerations  with  respect  to  time. 

The  basic  measuring  instrument  of  an  inertial  navigation 
system  is  the  accelerometer,  an  instrument  which  measures 
acceleration  along  a single  axis.  Inside  the  accelerometer 
is  a pendulous  mass  which  is  free  to  rotate  about  a pivot 
axis  in  the  instrument.  There  is  an  electric  pickoff  which 
converts  the  rotation  of  the  mass  about  the  pivot  axis  to  an 
output  signal.  The  output  signal  is  fed  to  a high  gain  amp1 f - 
fier  and  the  output  of  the  amplifier  is  connected  back  to 
torquing  coil  on  the  accelerometer.  When  an  acceleratioi 
present,  a current  is  sent  back  to  the  torquer  which  is  p 
cisely  the  amount  required  to  restore  the  mass  to  its  ini 
position.  The  accelerometer  torquer  can  be  restored  witl 
direct  current  or  it  can  be  restored  with  pulses.  Most 
modern  systems  use  pulse  restoration  because  of  the  ease 
which  the  accelerometer  output  can  be  processed  by  a dig^ 
computer.  A simplified  block  diagram  is  shown  in  Figure 
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The  current  fed  to  the  torquer  is  proportional  to  the 
measured  acceleration  and  provides  the  electrical  signal 
which  is  fed  to  the  navigation  computer.  The  computer 
integrates  the  acceleration  to  produce  velocity  and  then 
integrates  the  velocity  to  compute  distance.  If  two 


Defined  in  glossary 
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accelerometers  are  mounted  at  right  angles  to  each  other 
on  a platform  which  is  maintained  level  with  respect  to 
the  earth  and  if  one  of  the  accelerometers  is  directed  to 
true  North,  it  is  possible  to  determine  the  distance  that 
the  platform  traveled  in  the  North-South  and  East-West 
directions.  If  the  platform  is  "told"  where  it  is  initial- 
ly, it  can  determine  where  it  is  on  the  face  of  the  earth 
at  all  times.  A simple  block  diagram  is  shown  in  Figure  2. 


Figure  2 AcceUremetur*  on  Plotfonn 


The  accelerometers  must  be  maintained  level  and  in 
proper  azimuth^or ientation  regardless  of  aircraft  attitude 
because  even  a slight  misalignment  can  introduce  serious 
errors.  This  proper  attitude  and  azimuth  orientation  is 
maintained  by  mounting  the  accelerometers  on  a gimbaled 
platform  with  gyroscopes  used  as  sensing  elements  to  con- 
trol platform  orientation.  A platform  which  is  controlled 
by  gyros  in  this  manner  is  referred  to  as  a stable  element. 

A stable  element  requires  three  gyros  mounted  at  right 
angles  to  each  other.  One  gyro  senses  movement  about  pitch, 
one  about  roll,  and  one  about  azimuth.  Some  platforms  only 
require  two  gyros  but  they  are  two  degree  of-  freedom  types; 
i.e.,  they  are  sensitive  about  two  axes.  The  stable  element 
is  mounted  on  gimbals  to  isolate  it  from  angular  motions  of 
the  aircraft.  The  operation  of  the  gimbal  driving  system  is 
illustrated « i«  Figure  3. 
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Figure  3 Simplified  Single-Axle,  Gyre-Stabilized  Platform 


The  high  speed  spinning  wheel  or  rotor  in  Figure  3 resist 
any  motion  and  exhibits  an  important  gyro  characteristic  call 
ed  rigidity.  Any  movement  about  the  bearing  axis  in  Figure  3 
generates  a pickoff  signal  which  is  amplified  and  drives  the 
stable  element  in  an  equal  and  opposite  direction  to  the  move 
ment . 

A practical  inertial  system  requires  that  the  platform  be 
stabilized  in  all  three  axis  of  operation  in  order  to  retain 
its  level  orientation  regardless  of  the  maneuvers  made  by  the 
aircraft.  Figure  4 illustrates  a four-gimbal  platform  config 
ration  as  actually  used  in  an  inertial  system.  The  extra  rol 
gimbal  is  provided  to  prevent  the  occurrence  of  a condition 
known  as  gimbal  lock  during  certain  aircraft  maneuvers.  The 
gimbals  are  oriented  so  that  aircraft  attitude  and  heading^- 
may  be  sensed  by  measuring  angles  between  the  gimbals  and  the 
platform  frame.  Synchros  transmit  this  information  to  the 
attitude  indicator  and  other  systems  in  the  aircraft. 
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The  inertial  navigation  system  described  up  to  this 
point  is  only  capable  of  navigating  on  a flat  non-moving 
earth.  This  is  because  the  gyros  are  stablized  with 
respect  to  inertial  space  and  therefore  will  cause  the 
stable  element  to  rotate  with  respect  to  the  earth  as 
the  earth  rotates  on  its  axis.  This  is  undesirable  for 
navigation  because  the  accelerometers  will  not  remain  level 
with  respect  to  the  direction  of  gravity.  This  characteris 
tic  is  illustrated  in  Figure  5. 


Gyro  wheel 


PLATFORM 


SOUTH 
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When  that  stable  element  is  mounted  in  an  aircraft  the 
same  phenomenon  occurs  at  the  combined  rates  of  the  earth 
rotation  plus  the  aircraft  velocity.  Xf  the  aircraft  flies 
North  so  as  to  remove  the  rotation  of  the  earth  from  the 
sensitive  axis  of  the  stable  element,  the  aircraft  "sees” 
a continuing  pitch  maneuver.  At  the  pole,  instead  of  the 
platform  being  level  with  the  surface  of  the  earth,  it  would 
now  be  tilted  90  degrees  off  level.  This  characteristic  is 
illustrated  in  Figure  6. 


NORTH  POLE 
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Obviously  this  is  not  a practical  system  of  navigation 
over  a spherical,  rotating  earth  so  some  method  of  keeping 
the  stable  element  perpendicular  to  the  earth's  gravity 
vector  is  required.  By  applying  a torque  to  the  appropriate 
axis  it  is  possible  to  maintain  the  stable  element  with 
respect  to  the  earth  and  aligned  in  asimuth  to  a known 
reference.  Gyro  torque  is  produced  by  sending  a current 
through  torquer  coils  attached  to  the  gyro  gimbal.  This 
torque  causes  the  gyro  to  process  at  a right  angle  to  the 
applied  torque.  The  computer  supplies  the  current  to 


properly  torque  the  gyros.  The  operation  of  the  platform 
with  proper  earth  rate  and  vehicle  rate  gyro  torquing  is 
illustrated  in  Figure  7. 
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The  torquing  necessary  to  maintain  the  stable  element 
level  with  respect  to  the  gravity  vector  is  referred  to  as 
the  transport  rate  torque.  The  mechanization  of  this  trans- 
port rate  requires  a computing  loop  which  is  said  to  be 
Schuler  tuned.  This  computing  loop  (which  is  the  same  for 
each  platform  axis,  except  azimuth)  is  shown  in  Figure  8. 
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The  torquing  rate  necessary  to  maintain  the  platform 
level  with  respect  to  the  surface  of  the  earth  is  equal  to 
liicrtf t velocity  divided  by  the  radius  of  the  earth.  The 
8yatem  is  mechanized  so  that  velocity  is  multiplied  by  the 
reciprocal  of  the  earth's  radius  as  shown  in  Figure  9. 


Ffeni  ScM«-T«tag 

When  an  inertial  platform  is  mechanized  in  this  way, 
its  response  is  similar  to  a pendulum  with  a length  of  3,400 
miles  and  a period  of  84.4  minutes  as  shown  in  Figure  10. 

In  effect,  the  center  of  gravity  of  the  pendulum  bob  re- 
mains at  the  center  of  the  earth  and  the  point  of  suspension 
is  at  the  aircraft.  The  point  of  suspension  of  a pendulum 
can  be  moved  without  causing  the  pendulum  to  oscillate  so 
Schuler  tuning  allows  the  platform  to  be  moved  about  the 
surface  of  the  earth  without  disturbing  the  platform  level. 


PENDULUM  OSCILLATION  PERIOD 
IF  DISTURBED  »|4.4  MINUTES 
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An  Important  advantage  of  Schuler  tuning  la  that  many  of 
the  Instrument  errors  are  constrained  from  lncreaalng  with 
time  and  Instead  are  oscillatory  In  their  buildup.  An 
example  of  a constrained  acceleration  error  Is  Illustrated 
In  Figure  11.  This  error  Is  due.  to  an  Initial  level  axis 
misalignment  of  approximately  one  minute  of  arc.  Note  that 
the  platform  la  Initially  tilted  to  the  left  so  the  accelero- 
meter sensea  a component  of  gravity.  The  system  "thinks"  It 
Is  moving  and  a velocity  Is  developed  which  torques  the  plat- 
form In  the  opposite  direction  which  tends  to  cancel  the 
original  tilt  error.  The  system  then  overshoots  level  and 
develops  an  acceleration  error  in  the  opposite  direction. 

This  oscillation  will  continue  until  the  system  Is  switched 
off  or  damped  by  a velocity  input  from  an  outside  source 
and  is  referred  to  as  a Schuler  oscillation. 
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Correcting  the  platform  for  transport  rate  error  allows 
the  system  to  be  used  for  navigation  but  not  very  accurately. 
There  are  three  more  sources  of  error  that  must  be  corrected 
before  the  system  will  navigate  accurately.  They  are 
the  Coriolis  effect,  the  oblateness  (flattening  at  the  poles) 
of  the  earth,  and  centripetal  acceleration. 


IT 


The  Coriolis  effect  exists  because  the  earth  rotates.  If 
the  earth  were  not  rotating,  a vehicle  flying  a straight 
ground  course  from  the  equator  to  the  North  pole  would  make 
a straight  track  as  seen  In  Figure  12a.  However,  because 
the  earth  is  rotating,  an  observer  in  space  looking  down  on 
the  earth  would  see  that  the  airplane  really  has  to  fly  a 
curved  track  In  space  In  order  to  make  the  desired  straight 
ground  track  over  the  earth  (refer  to  Figure  12b).  Relative 
to  space,  the  airplane  must  continuously  change  the  magni- 
tude and  direction  of  its  tangential  velocity. 

II  n = EARTH  ROTATION  RATE  - 15.04  4e«/lir 
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Regardless  of  the  direction  of  the  vehicle's  horizontal 
velocity,  the  Coriolis  effect  appears  to  be  an  acceleration 
to  the  right  In  the  Northern  Hemisphere  and  to  the  left  in 
the  Southern.  Although  the  magnitude  of  the  Coriolis  accele- 
ration is  small,  it  can  cause  a significant  navigation  error 
when  flight  times  are  long.  To  compensate  for  the  Coriolis 
effect,  a correction  factor  is  introduced  into  the  output 
signals  of  the  accelerometers  which  is  equal  to  2ftvsln$, 
where  0 is  equal  to  earth  rate  (15.04  deg/hr),  v is  velocity 
(N-S  or  E-W)  , and  ♦ is  latitude.^* 
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The  oblateness  of  the  earth  produces  a spurious  accelera- 
tion In  the  N-S  accelerometer  when  the  plumb  line  to  the 
center  of  the  earth  does  not  exactly  coincide  with  the  true, 
vertical  as  shown  in  Figure  13.  The  navigation  computer  must 
provide  an  earth  radius  correction  term  to  the  N-S  accelero- 
meter to  correct  for  the  oblateness. 


Centripetal  acceleration  (ac)  is  a true  acceleration  with 
respect  to  inertial  space  that  is  created  by  a vehicle  travel- 
ing about  a spheroid  on  any  course  except  a great  circle 
course.  Figure  14  shows  an  aircraft  flying  East  along  a line 
of  constant  latitude  of  radius  r'.  The  center  of  curvature 
of  the  path  of  the  aircraft  does  not  pass  through  the  center 
of  the  earth  so  a centripetal  acceleration  is  generated.  Since 
the  path  shown  is  only  to  the  East  and  is  in  the  Northern 
hemisphere,  a South  acceleration  component  is  generated. 

To  correct  for  the  centripetal  acceleration,  a correction 
factor  equal  to  -v2/r'tan*  must  be  added  to  the  output 
signal  of  the  North-South  accelerometer. 


So  far  in  this  discussion,  only  the  local  vertical.  North 
pointing  inertial  navigation  system  has  been  considered. 

This  type  of  INS  is  referred  to  as  a semi-analytical  inertial 
system  and  is  the  most  common  system  in  use  today.  The 
platform  gimbal  structure  is  simple  and  the  computer  mech- 
anization is  easy.  The  semi-analytical  system  maintains 
the  stable  element  perpendicular  to  the  earth's  gravity 
vector  at  all  times.  Accelerometer  outputs  are  con- 
verted to  velocity  and  distance.  Velocities  are  used  to 
torque  the  stable  element  to  maintain  the  platform  normal 
to  the  earth  reference.  The  stable  element  is  aligned 
in  azimuth  to  an  azimuth  reference.  A simple  block  diagram 
of  a semi-analytical  inertial  system  is  shown  in  Figure  15. 
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One  of  the  major  disadvantages  arid  sources  of  navigation 
error  of  the  semi-analy tical  Inertial  platform  is  the  re- 
quirement to  torque  the  gyros  to  maintain  the  platform  per- 
pendicular to  the  earthJ  s gravity  vector.  This  problem  can 
be  overcome  by  using  an  analytical  navigation  system  which 
does  not  require  gyro  torquing.  Because  the  analytical 
navigation  platform  remains  fixed  in  space  and  rotates  with 
respect  to  earth,  the  accelerometers  must  sense  the  gravity 
component  and  vehicle  accelerations.  For  navigation  purposes 
only  the  vehicle  accelerations  are  used  and  the  gravitational 
accelerations  must  be  cancelled  out.  Calculating  the  earth's 
gravitational  acceleration  is  an  extremely  difficult  problem 
and  requires  that  an  enormous  amount  of  data  be  stored  in 
the  computer  memory.  Acceleration  corrections  and  present 
position  information  are  developed  in  the  computer  and  coordi- 
nate converter.  A simple  block  diagram  of  the  analytical 
inertial  system  is  shown  in  Figure  16. 
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A third  type  of  inertial  navigation  system  is  the  strap- 
down  system.  As  the  name  implies,  the  strap-down  system  does 
not  have  a conventional  gimbal  mounted  stable  element.  In- 
stead, the  gyros  and  accelerometers  are  mounted  directly  on 
the  vehicle  frame.  The  computer  in  the  strap-down  system 
must  compute  the  B matrix  computations  necessary  to  specify 
vehicular  attitude  with  respect  to  an  inertial  reference 
frame.  This  computation  is  usually  in  direction  cosine 
notation,  direction  cosines  being  any  space  vector  represent- 
ed by  three  cosines. 

The  coordinate  converter  utilizes  inputs  from  the 
accelerometers  and  the  B matrix  to  resolve  accelerations 
in  an  inertial  reference.  A position  computer  converts 
inertial  accelerations  and  altitude  information  to 
cartesian  coordinates  representing  vehicle  position  in 
Inertial  space.  A vector  computation  then  provides 
outputs  in  latitude  and  longitude.1 


Severe  torquing  requirements  are  placed  on  the 
gyros  in  a strap-down  system  because  the  gyros  change 
attitude  at  the  same  rate  as  the  vehicle.  A simple  block 
diagram  of  a strap-down  inertial  system  is  shown  in  Figure 
17. 
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ERROR  ANALYSIS 


GENERAL 

i 

The  magnitude  of  the  errors  in  an  inertial  system  may 
be  considered  as  random  variables  at  any  specific  flight 
time.  In  the  analysis  that  follows,  it  is  assumed  that  all 
errors  are  independent,  linear,  and  with  negligible  latitude 
effects  on  Coriolis  and  earth  rate  terms. 

ERROR  SOORCES 

Gyro  : 

Bias  (drift) 

Proportional  bias  (g  sensitive)1 
Anisoelasticity  (g2  sensitive)1 
Torquing  error 
Random  drift 
Scale  factor 
Accelerometer : 


Bias 

Scale  factor 
Cross  coupling1 
Vibropendulous 1 
Platform  : 

Initial  level 

Initial  azimuth  alignment 

Servo  errors 

Component  non-orthogonality 
Gyro 

Accelerometer 
defined  in  glossary 
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Computer 


Round  off  error  (register) 

Truncation  error  (integration  approximation) 
Commutation  error  (direction  cosine) 

Pick  off  error  (acceleration  roundoff  error) 

General : 

Geophysical  data 
Guidance  equations 
Target  locations 

The  local  vertical  system  will  be  studied  in  detail  be- 
cause it  is  a conventional  system  that  is  relatively  linear 
and,  consequently,  can  be  handled  with  Laplace  transforma- 
tion. The  understanding  of  the  local  vertical  system  is  basic 
to  all  systems  operating  near  the  earth.  As  shown  in  Figure 
18  and  19,  errors  can  originate  from  various  sources.  For 
this  discussion  only  error  sources  of  significant  magnitude 
(greater  than  2,000  feet/hour)  will  be  considered  and  all  error 
sources  will  be  considered  to  be  step  functions. 
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Figure  19  Errer  Sources  in  Azimuth  Axis 


VELOCITY  ERRORS 

A velocity  error  (ev)  will  cause  the  platform  to  torque 
out  of  level.  The  out  of  level  platform  will  sense  a com- 
ponent of  gravity  which  it  interprets  as  an  acceleration.  This 
signal  is  integrated  into  a velocity  of  opposite  polarity  to 
the  initial  velocity  error.  These  velocities  oscillate  at  the 
characteristic  Schuler  period  of  84  minutes.  Refer  to  Figure 
18  where  ev  is  a velocity  error  inside  the  Schuler  loop.  From 
Figure  18,  the  block  diagram  shown  in  Figure  20  can  be  con- 
structed . 


Figure  20  Velocity  Error  Loop 


The  transfer  function  of  this  error  in  the  system  is: 


ed  = S 
ev  S2+X2 


converting  to  time  domain 
sinXt 

E^  = ev  j Equation  1 

ev  = velocity  error  in  ft/sec 
e<j  ■ distance  error  in  ft 
g = gravity  constant  = 32.2  ft/sec2 
r ■ earth  radius  = 2.09xl07  ft 
t ■ time  in  seconds 

Equation  1 can  be  converted  to  a more  useful  form. 

= .13Evsinu)t  Equation  2 

Ed  = distance  error  in  miles 
a)  = 4.46  radians/hr 
t = time  in  hours 

The  velocity  errors  inside  the  Schuler  loop  generate 
position  errors  like  the  ones  shown  in  Figure  21.  Velocity 
errors  are  easily  diagnosed  if  the  inertial  velocities  can 
be  compared  with  actual  ground  reference  velocities. 
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ACCELERATION  ERRORS 


Acceleration  errors  (ea)  are  integrated  into  an  erroneous 
velocity  which,  thru  the  Schuler  loop,  torques  the  platform 
out  of  level.  The  out  of  level  accelerometers  sense  a com- 
ponent of  gravity  which  is  opposite  in  polarity  to  the  accel- 
eration error.  For  ease  of  analysis,  the  input  error  is  assumed 
to  be  a step  function.  The  vector  sum  of  the  acceleration 
error  and  the  gravity  error  oscillates  at  the  Schuler  frequency. 
The  error  in  computed  distance  is  the  double  integral  of  this 
acceleration  error.  In  block  diagram  form,  this  dynamic  con- 
dition is  shown  in  Figure  22. 


The  transfer  function  of  this  error  is 


Multiplying  by 


for  a step  input  gives 


Converting  to  the  time  domain 


Equation  3 


Equation  3 can  ba  converted  to  a more  useful  form. 


£4  ■ 3 .45e4(l-cosit>t)  Equation  4 

ea  ■ distance  error  in  miles 
ca  “ acceleration  error  in  ft/sec2 
w = 4.46  radians/hour 
t s time  in  hours 

Acceleration  errors  generate  position  errors  as  shown 
in  Figure  23.  Note  that  position  error  does  not  increase 
with  time  but  is  constrained  by  the  mechanization  of  the 
Schuler  loop.  Acceleration  errors  are  difficult  to  isolate 
in  any  kind  of  operating  environment  other  than  a laboratory 
because  the  accelerometer  outputs  are  rarely  available 
external  to  the  inertial  platform  for  measurement  by  the 
instrumentation  system.  The  normal  method  of  examining  the 
accelerometer  outputs  is  to  differentiate  the  velocity 
outputs. 
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LEVEL  GYRO  DRIFT  ERRORS 

Level  gyro  drift  (e&)  is  one  of  the  most  common  sources 
of  inertial  error.  The  error  is  usually  the  result  of  an 
improper  gyro  bias  being  applied  to  one  or  both  of  the  level 
axis  gyros.  The  block  diagram  for  the  error  is  shown  in 
Figure  24. 
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When  an  error  enters  the  system  as  a level  gyro  drift 
error  the  platform  is  unleveled.  The  accelerometers  pick 
up  a component  of  gravity  which  is  integrated  into  velocity. 
Through  the  Schuler  loop  the  erroneous  velocity  computation 
torques  the  level  stabilization  gyro  in  opposition  to  the 
error  signal.  However,  in  an  undamped  pure  inertial  system, 
the  velocity  computation  overshoots  the  error  signal  and 
oscillates  as  a negative  sine  wave  at  the  Schuler  frequency 
about  a ramp  function  of  time.  The  transfer  of  the  error  is: 


ed 

— ■ r 


S { Sz  + X 2 j 


Multiplying  by  1/S  for  a step  function  input  error  gives 


Ed  - e*r 


S2 (S2+X2) 





Converted  to  the  time  domain 


sinXt. 

Ed  “ etr(t ) 


Equation  5 


e<j  » distance  error  in  feet 
e*  * gyro  drift  in  radians/sec 
r - earth  radius  » 2.09  x 107  feet 
t = time  in  seconds 

Equation  5 can  be  converted  to  a more  useful  form  by  con- 
verting units. 

ea  * 60 e£  ( t-0 . 22s inoit ) Equation  6 

* gyro  drift  in  deg/hr 
t * time  in  hours 
at  * 4 .46  radians/hr 
ed  = distance  error  in  miles 

Level  axis  gyro  drift  generates  position  errors  as  shown 
in  Figure  25. 


AZIMUTH  DRIFT  ERRORS 


The  most  significant  cross-coupling  error  between  axis 
is  that  of  azimuth  drift  rate  (6a) . This  error  (Figure  19) 
is  integrated  by  the  azimuth  stabilization  loop  resulting 
in  an  azimuth  misalignment  angled  When  misaligned,  the  east 
gyro  picks  up  a component  of  the  earth's  rate  of  rotation 
which  torques  the  stable  element  out  of  level.  The  re- 
sulting gravity  error  sensed  by  the  accelerometer  is 
integrated  into  velocity  and  torques  the  platform  out  of 
level  opposite  to  the  earth  rate  torque  through  the  Schuler 
loop.  The  computed  velocity  signal  overshoots  and  oscil- 
lates about  the  earth's  rate  torque.  The  block  diagram 
for  this  error  is  shown  in  Figure  26. 
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From  this,  the  error  transfer  function  is  derived. 
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Multiplying  by  1/S  for  a step  function  arror  input  gives 


*4  * «.‘>co»«r(s5(sJ*'J)J 
Converting  to  the  time  domain 

ed  - «aftcos*r(f  - Equation  7 

ed  ■ distance  error  in  feet 
6a  » azimuth  drift  rate  in  radians/sec 
ft  - earth  rotation  rate  - 15.04  deg/hr 
4 - local  latitude 
r « earth  radius-2.09  x 107  ft 
t - time  in  seconds 

Equation  7 can  be  reduced  to  a more  useful  form. 

ed  • 7.86cos#fiat2  Equation  8 

ed  ■ distance  error  in  miles 
6a--  azimuth  gyro  drift  in  deg/hr 
* - local  latitude 
t - time  in  hours 
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This  error  is  predominately  in  the  North-South  direction 
because  the  effect  of  earth  rate  upon  the  north  gyro  will 
be  in  error  by  the  cosine  of  the  aximuth  misalignment.  The 
cosine  does  not  change  significantly  for  small  error  angles. 
The  distance  resulting  from  aximuth  gyro  drift  is  shown  in 
Figure  27.  Note  that  the  error  is  in  the  North-South  axis 
only  and  that  it  increases  in  magnitude  as  the  square  of  time. 
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VELOCITY  ERRORS  OUTSIDE  THE  SCHULER  LOOP 

Velocity  errors  entering  from  outside  the  Schuler  loop 
(evo  in  Figure  18)  affect  the  rate  of  the  second  integration 
and  cause  the  distance  computation  error  eg  to  have  a straight 
line  increase  when  the  error  is  constant. 


INITIAL  AZIMUTH  MISALIGNMENT 


Initial  aximuth  misalignment  (e0  in  Figure 
to  error  in  resolving  horizonal  accelerations, 
duces  errors  by  causing  the  east  gyro  to  sense 


19)  contributes 
It  also  intro- 
earth  rate. 


INITIAL  LEVEL  MISALIGNMENT 

In  addition  to  accelerometer  null  uncertainties,  there 
may  be  acceleration  errors  because  of  initial  level  misalign- 
ment errors  (e©,  e<j>  in  Figure  18)  due  to  a dead  band  in  the 
servo  null  of  the  platform  gimbals.  This  error  may  be  written 
as  ca  - gsince  for  a level  misalignment  in  the  pitch  axis 
where  g is  the  acceleration  due  to  gravity  in  f eet-per-second- 
per-second  and  e@  is  the  level  misalignment  of  the  pitch  axis 
in  radians.  An  error  in  the  roll  axis  would  be  written  as 
ea  - gsine^ . 
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TEST  PLANNING 


ADVANCE  PLANNING 

. When  the  Flight  Test  Center  is  assigned  as  the  Responsible 

Test  Organization  for  any  flight  testing,  a certain  amount  of 
advance  planning  is  required.  If  the  testing  is  to  include 
an  inertial  system,  an  estimate  of  the  number  of  flights  and 
flying  hours  required  to  evaluate  the  system  will  be  made. 

The  evaluation  of  the  INS  should  be  divided  into  two 
separate  areas,  accuracy  and  operational  suitability.  To 
predict  the  accuracy  of  an  INS  with  an  85  percent  confidence 
in  the  validity  of  that  prediction  requires  a minimum  of  eight 
dedicated  and  valid  flights  on  an  Instrumented  aircraft  for 
each  alignment  mode  that  has  an  accuracy  specification  written 
against  it.  The  eight  flights  per  alignment  mode  produce  a 
Circular  Error  Probable^(CEP)  prediction.  Each  flight  length 
should  be  in  excess  of  one  Schuler  period  of  one  hour  and  24 
minutes:.  It  is  desirable  to  exceed  two  Schuler  periods  but 
that  is  not  a requirement  as  long  as  the  test  flight  lengths 
are  compatible  with  the  operational  mission  requirements  of 
the  aircraft.  Operational  suitability  does  not  require  dedi- 
cated test  flights  but  adds  0.1  flying  hours  to  every  test 
flight. 

. 3 

To  determine  the  Inflight  accuracy  of  an  inertial  system 
with  a reasonable  confidence  level  requires  an  on-board  instru- 
mentation system.  The  inertial  parameters  must  be  measured 
with  a degree  of  precision  that  will  allow  testing  to  the 
defined  requirements.  The  specific  parameters  to  be  measured 
for  each  inertial  system  will  change  somewhat  from  aircraft 
to  aircraft  but  the  general  requirements  do  not  change.  The 
following  list  is  the  general  measurements  required  for  quanti- 
tatively testing  an  inertial  system: 

Inertial  longitude 
Inertial  latitude 
Inertial  elevation 
Inertial  ground  speed 
Inertial  ground  track 

/ Inertial  velocity  in  N-S  direction 

Inertial  velocity  in  E-W  direction 
Inertial  vertical  velocity 
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Inertial  wander  angle 
Inertial  roll  angle 
Inertial  pitch  angle 

Inertially  computed  magnetic  variation 
Inertially  computed  wind  speed 
Inertially  computed  wind  direction 
True  airspeed  input  to  inertial  system 
Magnetic^ heading  input  to  inertial  system 
Pressure  altitude  input  to  inertial  system 
Roll  angle  from  auxiliary  reference  system 
Pitch  angle  from  auxiliary  reference  system 
Magnetic  heading  from  auxiliary  reference  system 
Radar  range 
Radar  bearing 
TACAN  range 
TACAN  bearing 

Present  position  update  command 
Enter  visual  fix  command 
Aircraft  angle  of  attack 

Aircraft  vertical  acceleration  at  the  center  of  gravity 
Tone  and  event  command 


As  soon  as  the  location  of  the  test  facilities  are 
defined  at  Edwards  AFB , the  systems  test  engineer  should 
determine  if  the  alignment  coordinates  and  local  magnetic 
variation  is  available  for  all  of  the  proposed  parking  spots 
for  the  test  aircraft.  The  local  magnetic  variation  should 
be  remeasured  for  each  of  the  packing  spots  approximately 
every  four  years.  This  is  because  the  magnetic  variation 
in  the  Edwards  AFB  area  changes  at  the  rate  of  1.5  minutes 
per  year  or  0.1  degrees  every  four  years. 

DETAILED  PLANNING 


Test  Information  Sheet  (TIS? : 


After  the  preliminary  test  planning  has  been  calcu- 
lated, documented,  and  provided  to  the  project  engineer,  the 
systems  analyst  (test  engineer  responsible  for  the  analysis  of 
the  flight  test  data  on  the  inertial  system)  should  begin  the 
detailed  test  planning.  This  planning  is  documented  on  a 
for  inclusion  in  the  test  plan.  The  TIS  form  shown  in  Figure 
28  has  the  major  topics  identified  that  are  to  be  included  in 
the  documentation.  These  topics  will  be  discussed  individually 
in  the  following  paragraphs: 

^Defined  in  glossary 


*•*  \ 


TITH  OF  TOT 
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Tt*  TVFt  LOCATION  OF  TOT 

□ fla*  □FNocaouNAL  Edvards  AFB 


T0TINS  ACTIVITY 

AFTTC 


INASANOOUI/UNUfUAL  TOT 

No 


1.0  References:  The  first  reference  to  be  listed 

is  the  Air  Force  management  document  that  is  used  to  generate 
the  aircraft  or  system  contract.  This  document  will  detail 
the  operational  requirements  that  were  intended  to  be  satisfied. 
The  next  reference  document  will  be  the  contractors  system 
specification  document.  In  this  document,  the  contractor 
tries  to  quantify  the  operational  requirements  contracted 
for  by  the  Air  Force.  The  third  reference  document 
should  be  the  detailed  INS  specification  document  published 
by  the  manufacturer  of  the  INS.  The  fourth  reference  docu- 
ment should  be  the  aircraft  flight  manual  showing  how  to 
operate  the  INS.  All  of  the  maintenance  documents  used  at 
the  AFFTC  to  maintain  the  system  should  also  be  referenced. 

2.0  Test  Item  Description:  The  description  should 

not  be  of  an  INS  in  general  but  rather,  what  makes  up  this 
particular  inertial  system.  Identify  all  of  the  aircraft 
components  that  are  considered  to  be  a part  of  the  inertial 
system  and  will  be  included  in  the  evaluation.  Describe  the 
alignment  modes  and  operating  modes.  Briefly  discuss  any 
software  that  will  be  considered  a part  of  the  INS. 

3.0  Test  Objective:  The  test  objective  will  be 

much  the  same  for  any  INS  evaluated.  It  will  normally 

be  "To  determine  the  operational  accuracy  and  usability  of 
the  Inertial  Navigation  System." 

4.0  Success  Criteria:  The  accuracy  portion  of  the 

testing  will  be  complete  when  sufficient  test  points  are 
available  to  give  a reasonable  estimate  of  the  INS  accuracy. 

This  estimate  requires  a minimum  of  eight  valid  data  flights 
for  each  mode  of  operation  and/or  alignment.  If  any 
development  changes  are  made  to  the  hardware  or  software 
after  the  accuracy  testing  has  been  accomplished  and  if  those 
changes  could  affect  system  accuracy,  then  sufficient  re- 
testing must  be  accomplished  to  demonstrate  the  impact  of 
the  change.  Operational  suitability  testing  will  continue 

on  a ride-along  basis  until  the  end  of  the  test  program. 

5.0  Data  Requirements:  The  end  point  data  for  all 

valid  data  flights  is  normally  presented  on  a circular  error 
plot  with  the  vertical  axis  showing  North-South  error  and 
the  horizontal  axis  showing  East-West  error.  The  accuracy 
requirement  is  shown  as  a circle  whose  radius  is  equal  to 
the  specification  CEP.  Each  flight  terminal  error  point  is 
normalized  to  the  time  duration  of  the  specification  and 
plotted.  The  average  value  of  all  the  plotted  points  is 
calculated  and  plotted.  The  predicted  CEP  based  on  the  test 
data  is  derived  and  plotted  according  to  a formula  that  is 


used  by  the  AFFTC  to  evaluate  all  inertial  systems  tested 
at  the  Center.  If  the  contractor  uses  a different  technique 
for  calculating  CEP,  that  value  should  also  be  calculated  and 
plotted . 

If  in-flight  accuracy  data  are  obtained,  the  data  will 
be  presented  as  an  error  time-history  plot  for  each  flight. 
Multiple  flights  in  an  alignment  mode  will  be  used  to  generate 
a CEP  time-history  plot.  A separate  CEP  plot  will  be  made  for 
each  alignment  mode  evaluated.  An  85  percent  confidence  limit 
should  be  shown  on  each  CEP  time-history  plot  along  with  the 
specification  for  the  INS. 

Ground  test  data  should  be  obtained  and  compared  with  the 
specifications  called  out  in  the  maintenance  technical  publi- 
cations for  ground  checking  the  platform.  The  relationship 
of  ground  accuracy  and  in-flight  accuracy  should  be  docu- 
mented . 


6.0  Test  Procedures:  This  section  will  contain  the 

overall  test  plans  for  the  inertial  platform.  It  should  be 
as  detailed  as  possible  at  the  time  the  TIS  is  prepared  and 
should  start  with  ground  testing  and  conclude  with  in-flight 
testing.  Specific  test  techniques  will  be  discussed  for 
various  types  of  testing  later  in  this  report. 

7.0  Support  Requirements:  This  section  should  detail 

all  of  the  services  furnished  by  organizations  other  than 
Systems  Engineering.  included  should  be  a list  of  the  instru- 
mentation parameters,  range  support,  data  support,  radar  track- 
ing support,  and  photographic  support. 

Data  Collection: 

Collection  of  the  INS  data  is  the  responsibility  of 
the  systems  test  engineer  assigned  to  evaluate  the  INS  system. 
The  test  engineer  must  plan  in  great  detail  how  the  test 
points  are  to  be  flown.  After  he  has  planned  how  to  obtain 
the  test  points  on  each  flight,  he  should  make  up  the  test 
cards  for  this  part  of  the  testing.  He  should  discuss 
the  test  cards  with  the  Project  Engineer  to  make  sure  that 
all  test  points  and  cards  are  practical.  A test  support 
summary  should  be  written  up  for  each  test  flight  detailing 
flight  speed  and  length,  photo  and/or  safety  chase  require- 
ments, radar  tracking  requirements,  photo  theodolite  tracking 
requirements  and  instrumentation  support  requirements.  This 
planning  must  be  complete  before  scheduling  a test  mission 
(approximately  two  weeks  before  flight) . 
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Instrumentation  support  requirements  planning  should 
produce  a prioritized  list  of  parameters  which  must  be  opera- 
ting for  data  analysis  on  each  flight.  The  list  should  in- 
clude measurements  of  aircraft  configurations,  aircraft 
attitude,  and  aircraft  dynamic  conditions  so  that  it  can  be 
determined  what  the  dynamic  conditions  of  flight  were  at 
specific  times  during  the  flight.  The  planning  should  also 
identify  any  special  instrumentation  pre-flight  and  post 
flight  requirements. 

A complete  history  will  be  kept  on  each  inertial 
platform  in  the  test  program.  When  an  inertial  platform 
is  initially  assigned  to  the  test  program,  a platform  history 
log  like  the  one  shown  in  Figure  29  is  initiated.  Every 
hour  of  platform  operation  during  the  test  program  should 
be  documented,  regardless  of  the  reason  for  the  operation. 

This  requires  the  cooperation  of  the  ground  maintenance 
personnel.  Control  procedures  must  be  worked  out  with  the 
maintenance  personnel  and  the  project  engineer  on  each 
test  program  to  make  sure  that  none  of  the  data  is  lost. 

The  accuracy  of  the  INS  platform  boresight  should  be 
determined  and  documented  at  some  time  during  the  test  program. 
This  is  necessary  to  prove  that  INS  data  obtained  on  any  parti- 
cular aircraft  is  valid  and  does  represent  what  can  be  expected 
out  of  the  using  command  aircraft  providing  they  are  also  within 
tolerance  on  the  boresight  alignment.  The  boresight  check  should 
be  made  as  early  as  possible  in  the  test  cycle.  If  the  plat- 
form is  out  of  tolerance,  all  data  obtained  prior  to  correcting 
the  boresight  is  invalid  and  can  not  be  used  to  demonstrate 
the  system  capability. 

Inertial  data  are  obtained  during  all-weather  testing 
and  should  be  analyzed.  However,  that  test  requirement  should 
be  considered  as  being  in  addition  to  the  baseline  data  that 
will  be  obtained  at  Edwards  AFB.  The  reason  for  this  is  that 
if  there  are  any  observable  differences  between  the  data  be- 
cause of  the  difference  in  environment,  there  must  be  enough 
data  at  each  condition  to  prove  and  present  the  difference. 

Keep  in  mind  that  the  all-weather  test  aircraft  will  be  tested 
at  a remote  site  so  it  may  be  impossible  to  analyze  the  data 
in  a very  timely  manner.  In  fact,  there  is  a good  probability 
that  the  test  engineers  that  accompany  the  all-weather  teat 
aircraft  to  a remote  test  site  will  not  have  the  background 
to  analyze  the  INS  test  data.  Their  primary  responsibility 
is  to  conduct  the  test,  not  to  analyze  subsystem  performance. 
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Ground  Testing; 

Sufficient  ground  testing  should  be  accomplished  to 
establish  the  relationship  of  the  ground  accuracy  to  the 
in-flight  accuracy  and  to  validate  the  ground  maintenance 
procedures  called  out  in  the  T.O.’s.  This  testing  should  be 
accomplished  by  the  engineer  responsible  for  analyzing  the 
INS  performance  but  could  be  accomplished  by  others.  As  an 
example,  maintenance  will  often  perform  a ground  drift  run  to 
check  out  the  system.  The  engineer  should  obtain  and  use  the 
data  from  the  drift  run. 

The  length  of  a ground  drift  run  should  never  be 
less  than  84  minutes.  A minimum  of  eight  ground  drift  runs 
should  be  made  for  each  alignment  mode  during  the  course  of 
the  test  program  to  obtain  a confidence  level  of  at  least 
85  percent  in  the  validity  of  the  data.  If  in-flight 
accuracy  data  are  obtained,  it  is  a good  policy  to  obtain 
drift  run  data  before  and  after  the  flight  using  the  same 
alignment  mode  as  for  the  flight.  However,  the  platform 
should  be  allowed  to  cool  down  to  ambient  temperature 
between  shutdown  and  the  next  alignment  and  that  requires 
a minimum  of  12  hours.  It  is  also  preferable  to  start  the 
alignment  on  a different  heading  than  the  shutdown  heading. 
This  will  insure  that  the  platform  will  have  a reasonable 
amount  of  shutdown  error  to  correct  during  alignment.  Align- 
ing on  various  headings  helps  to  identify  the  effectiveness 
of  the  heading  bias  correction  constant.  A sample  INS 
ground  drift  run  card  is  shown  in  Figure  30. 


Operational  Testing: 


For  operational  testing,  the  data  that  will  be  plotted 
is  normalized  endpoint  data.  The  inertial  platform 
history  log  shown  in  Figure  29  has  all  of  the  data  required 
to  produce  endpoint  data  plots.  The  flight  data  on  that  log 
will  be  obtained  by  the  systems  test  engineer  from  the 
flight  crew  at  postflight  debrief.  Sample  alignment  and  shut- 
down cards  are  shown  in  Figures  31  and  32.  These  cards  should 
be  prepared  by  the  systems  test  engineer  and  given  to  the 
project  engineer  at  least  one  day  prior  to  flight.  The 
REMARKS  entry  shown  on  Figure  31  is  to  document  the  wind  con- 
ditions during  the  alignment,  any  observable  vibration  or 
movement  of  the  aircraft  during  the  alignment,  and  whether  the 
power  source  was  switched  from  ground  power  to  aircraft  power 
during  the  alignment. 
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INERTIAL  ALIGNMENT  1 1 INERTIAL  SHUTDOWN 


Figure  31 


Inertial  Alignment  Test  Card 


The  flight  data  will  be  presented  on  data  plots  that 
reflect  the  alignment  mode.  To  demonstrate  any  type  of  rapid 
alignment,  the  INS  should  be  allowed  to  soak  at  ambient  con- 
ditions for  a minimum  of  12  hours. before  starting  the  align- 
ment. If  the  ambient  soak  conditions  can't  be  met,  then  the 
platform  should  be  aligned  in  the  standard  alignment  mode 
(gyro-compass  for  a semi-analytical  system) . 

If  the  inertial  platform  has  to  be  shut  down  in 
flight  for  any  reason,  the  data  recorded  before  platform 
shutdown  can  not  be  used  for  endpoint  accuracy  data.  How- 
ever, if  the  platform  is  brought  back  up  with  an  in-flight 

alignment,  that  data  can  be  used  as  end-point  data  for  the 
in-flight  alignment  accuracy  plot  providing  it  is  possible 
to  define  the  accuracy  of  the  aircraft  position  versus  the 
entered  initialization  coordinates. 

Often  during  a test  program  it  will  be  necessary  to 
enter  a present  position  correction  to  the  INS  during  flight. 

For  a non-instrumented  aircraft,  that  almost  always  prevents 
the  test  engineer  from  using  the  data  from  that  flight  for  any- 
thing other  than  a qualitative  description  of  how  well  it 
worked.  This  is  because  it  is  necessary  to  know  where  the 

aircraft  was  immediately  before  update,  how  much  error  was  in 

the  INS  at  that  time,  and  how  much  of  the  error  was  removed 
by  the  update  in  order  to  use  the  data  from  the  flight  in 
any  kind  of  meaningful  analysis.  All  of  this  information 
must  be  time  correlated  to  the  alignment  and  shutdown  data. 

Many  inertial  systems  tested  in  recent  years  are 
bounded-error  systems.  This  means  that  position  and  velocity 
parameters  from  the  INS  are  mixed  with  data  from  other  sources 
by  a statistical  filter  to  arrive  at  a computer  estimate  of 
position  and  velocity.  This  computer  estimate  is  then  fed 
back  to  the  INS  to  try  to  keep  the  difference  between  them  as 
small  as  possible.  It  is  not  really  possible  to  isolate  the 
INS  accuracy  because  the  accuracy  of  the  other  systems  show 
up  in  the  results. 

Air  Force  accuracy  requirements  are  nearly  always 
specified  for  operational  conditions.  However,  the  con- 
tractor normally  demonstrates  specification  compliance  under 
a limited  set  of  carefully  controlled  conditions.  The 
initial  responsibility  for  quantitative  testing  under 
operational  conditions  falls  to  the  Flight  Test  Center  and 
the  systems  test  engineer  should  lay  out  his  testing  require- 
ments with  that  in  mind.  Alignments  should  be  made  in  all 
available  modes  and  headings.  Maneuvers  should  be  specified 
on  a significant  number  of  flights  that  will  stress  the  inertial 
platform  with  the  same  type  loads  expected  operationally. 

The  systems  test  engineer  must  personally  prepare  the 
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alignment  and  shutdown  cards  for  each  flight  and  review 
all  flight  cards  to  make  sure  that  the  flight  will  represent 
a fair  operational  challenge  to  the  INS. 
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Inflight  Accuracy  Testing; 

The  data  objective  of  inflight  accuracy  testing  is 
to  produce  a single  plot  that  compares  the  predicted  CEP 
based  on  actual  test  data  with  the  specification  CEP  for 
each  alignment  mode  evaluated.  It  is  desired  to  have  a 
minimum  of  eight  flights  of  valid  data  for  each  mode  to 
calculate  a CEP  with  a reasonable  confidence  level. 

For  inflight  accuracy  testing  of  current  inertial 
systems,  an  instrumentation  system  is  usually  required  and 
is  always  highly  desirable.  The  system  should  record  the 
on-board  position  coordinates  with  respect  to  a time  that 
is  corr elateable  to  ground  tracking  data.  As  inertial 
accuracy  flights  are  dedicated  flights,  the  instrumentation 
recorder  should  be  turned  on  and  off  at  the  desired  test 
point  intervals.  An  instrumentation  event  tone  should  be 
transmitted  briefly  during  each  recorder-on  cycle  so  that 
each  segment  of  recorded  data  can  be  individually  checked 
for  time  correlation. 

Range  tracking  accuracy  is  an  important  considera- 
tion for  inertial  accuracy  testing.  For  systems  that  have 
been  tested  at  Edwards  in  the  past,  the  accuracy  of  the  space 
positioning  radar  has  been  sufficient  if  the  radar  grazing 
angle  is  kept  a minimum  of  two  degrees  above  the  horizontal 
plane  of  the  radar  antenna.  The  azimuth  angle  accuracy  of 
the  radar  is  specified  as  0.2  mils  which  is  1.2  feet  per  mile 
of  range.  Figure  33  is  a plot  of  the  radar  azimuth  accuracy 
as  a function  of  range  from  the  tracking  antenna.  The  minimum 
MSL  altitude  of  the  aircraft  required  to  obtain  this  accuracy 
is  shown  at  25  mile  increments.  The  ranging  accuracy  for 
the  radar  is  specified  as  ±12  feet  at  all  ranges. 

The  flight  profile  for  an  accuracy  flight  can  be  very 
flexible  but  because  of  the  amount  of  data  required,  it  is 
recommended  that  a rectangular  pattern  be  flown  along  the 
lines  of  constant  latitude  and  longitude  with  the  length 
of  the  legs  being  as  long  as  possible  without  exceeding  the 
altitude  and  accuracy  restrictions.  High  rate  turns  and  pitch 
maneuvers  should  be  programmed  into  the  individual  legs  to 
task  the  inertial  system  as  much  as  possible.  The  systems 
test  engineer  should  monitor  the  test  flight  from  whatever 
location  that  is  available  to  him  that  provides  him  the 
greatest  insight  to  the  data  being  produced. 
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The  data  to  be  collected  during  the  alignment  before 
an  inertial  accuracy  flight  is  shown  in  Figure  34.  Observe 
that  the  data  card  requires  the  crewmember  to  duplicate  many 
of  the  parameters  being  recorded  by  the  instrumentation  system. 
This  is  so  that  if  the  instrumentation  system  should  fail, 
there  would  be  sufficient  data  available  to  keep  the  inertial 
Platform  History  Log  current  and  use  for  end-point  accuracy 
data.  Initial  alignment  headings  should  be  varied  so  that  a 
minimum  of  two  alignments  will  be  performed  in  each  compass 
quadrant . 

The  next  data  point  is  on  the  taxi-way  just  before 
the  takeoff  end  of  the  runway.  There  is  an  inertial  checkpoint 
marked  on  the  taxi-way  at  both  ends  and  the  center  of  the 
runway  at  Edwards  AFB . Photographs  of  the  checkpoints  at 
the  ends  of  the  runway  are  shown  in  Figures  35  and  36.  When 
the  nosewheel  is  over  the  INS  checkpoint,  the  pilot  should  stop 
the  aircraft,  place  an  event  mark  on  the  instrumentation  tape, 
and  write  down  the  time  and  coordinates.  A suggested  flight 
card  format  is  shown  in  Figure  37. 

During  the  flight,  test  points  should  be  recorded 
every  five  minutes.  Because  of  the  number  of  test  points  in- 
volved, the  number  of  test  cards  for  inflight  accuracy  data 
could  become  unwieldy  for  an  extended  mission.  It  is  important 
to  lay  out  the  test  in  such  a way  as  to  assure  being  able  to 
correlate  data.  A sample  test  card  is  shown  in  Figure  38. 

As  stated  earlier  in  this  report,  the  flight  should  last  a 
minimum  of  one  full  Schuler  cycle  of  84  minutes  so  there 
should  be  a minimum  of  17  in-flight  test  points  five  minutes 
apart . 


The  Inertial  Platform  History  Log  should  be  filled  out 
by  the  systems  test  engineer  from  data  obtained  in  the  post- 
flight debriefing  just  as  if  it  was  a non-instrumented  flight. 
The  fact  that  in-flight  accuracy  data  exists  for  a flight  should 
be  noted  in  the  Remarks  column. 

REVIEW  OF  EXISTING  DATA 

Normally,  the  Flight  Test  Center  will  not  test  an  INS 
that  does  not  have  a considerable  amount  of  test  data  already 
available  on  it.  That  test  data  may  be  from  the  INS  manufac- 
turers development  testing,  from  testing  conducted  by  the 
Central  Inertial  Guidance  Test  Facility  (CIGTF)  at  Holloman 
AFB,  NM,  from  the  avionics  development  and  integration  test 
program  of  the  airframe  contractor,  or  from  test  reports  on 
other  aircraft  with  the  same  inertial  system  installed.  The 
systems  test  engineer  should  attempt  to  obtain  as  much  of  this 
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INERTIAL  ALIGNMENT 


1.  Align  mode 


2.  Instrunentation  recorder  ON 


3.  Function  switch  to  align 


4.  Enter  coordinates  and  reset 


5.  Heat  light  out 


6.  Record  outside  air  temperature 


7.  Alignment  complete 


8.  System  heading 


9.  System  magnetic  variation 


MAG  VAR 


10.  Select  NAV  position 


11.  Instrumentation  recorder  OFF 


Figure  34  Alignment  Test  Card  - Instrumented  Aircraft 
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Runway  22  at  Edwards  AFB 


Runway  4 at  Edwards  AFB 


Figure  36  Inertial  Accuracy  Checkpoint  No.  2 


INERTIAL  TEST  REQUIREMENTS  AT  TAKEOFF  I I INERTIAL  TEST  CARD 


data  as  possibla  prior  to  ths  start  of  tha  tast  program  and 
bo  praparad  to  compare  it  with  Flight  Tost  Cantor  data  as 
that  bacomas  available. 

Xf  tha  data  from  all  sources  does  not  correlate,  the 
reason  must  be  determined.  If  the  data  does  correlate,  it 
might  be  possible  to  reduce  the  number  of  data  flights  re- 
quired at  the  Flight  Test  Center.  However,  this  is  a function 
of  the  confidence  that  the  systems  test  engineer  has  in  the 
data,  so  reducing  the  number  of  test  flights  should  not  be 
included  in  any  of  the  planning. 

REVIEW  OF  AGE  AND  TECHNICAL  PUBLICATIONS 

It  is  the  responsibility  of  the  systems  test  engineer  to 
review  all  applicable  technical  publications  and  evaluate 
the  contractor  furnished  support  hardware  to  be  used  by 
Air  Force  maintenance  personnel.  The  technical  publica- 
tions should  give  clear  and  accurate  directions  on  opera- 
tion, troubleshooting  and  maintenance  procedures  for  the 
INS.  The  test  equipment  should  provide  assistance  in  the 
trouble-shooting  and  repair  of  the  system. 

An  adequate  evaluation  of  the  AGE  and  technical  publi- 
cations requires  the  assistance  and  cooperation  of  the 
maintenance  personnel  assigned  to  the  test  program.  The 
systems  test  engineer  should  try  to  keep  the  maintenance 
personnel  fully  aware  of  what  he  is  trying  to  accomplish 
and  the  results  from  data  already  gathered.  This  includes 
showing  them  the  data,  showing  how  it  correlates,  and 
discretely  pointing  out  the  mistakes  that  maintenance 
personnel  may  make  that  can  cause  the  data  to  be  invalid. 

Any  recommended  changes  in  the  maintenance  procedures  should 
be  discussed  with  the  maintenance  personnel  before  submitting 
to  the  SPO. 

REPORTS 

The  systems  test  engineer  should  plan  for  the  reporting 
phase  of  the  INS  testing  just  as  carefully  as  he  does  the 
rest  of  the  test  program.  The  reports  must  present  all 
knowledge  obtained  as  a result  of  the  testing. 

The  first  type  of  reporting  that  may  be  required  of  the 
test  engineer  is  for  Air  Force  Preliminary  Evaluations.  This 
type  of  testing  is  normally  a quick  overview  where  the  entire 
aircraft  is  evaluated  on  a very  small  number  of  missions. 

If  an  INS  had  gross  problems  that  would  make  it  operationally 
unacceptable,  this  type  of  testing  should  identify  it.  How- 


ever,  the  likelihood  of  detecting  subtle  design  problems  or 
accuracy  problems  is  not  very  great.  The  data  should  be 
presented  as  a written  summary  of  what  was  observed  along 
with  conclusions  about  the  system  operation. 

During  the  conduct  of  a flight  test  program,  the  Deficiency 
Report  (DR)  is  used  to  document  and  track  deficiencies  in  the 
inertial  system.  The  specific  directions  for  using  the  DR 
form  can  be  found  in  AFFTCR  80-2.  Generally  speaking,  the 
DR  is  used  to  identify  suspected  problems  to  the  SPO  and  the 
contractor.  The  systems  test  engineer  will  have  the 
responsibility  of  keeping  track  of  the  DR's  on  the  INS 
and  whatever  corrective  actions  that  take  place.  A DR  form 
is  shown  in  Figure  39. 

The  progress  report  is  a periodic  management  reporting 
tool  to  keep  the  report  addressees  informed  on  the  progress  of 
the  testing  and  to  provide  an  insight  into  any  problems  that 
exist.  The  progress  report  should  summarise  the  INS  testing 
for  the  report  period,  summarize  the  results,  relate  that 
to  what  has  been  conducted  previously,  and  describe  what 
remains  to  be  accomplished.  A summary  of  the  deficiency 
reports  submitted  during  the  reporting  period  should  be  in- 
cluded. If  several  flights  are  involved,  an  INS  flight  log 
like  the  one  shown  in  Figure  40  should  be  included.  Note  that 
all  of  the  information  presented  on  the  flight  log  can  be 
obtained  from  the  Inertial  Platform  History  Log  (Figure  29) . 

The  final  report  should  cover  everything  that  is  known 
by  the  Test  Center  about  the  INS.  The  body  of  the  report 
should  present  summarized  data  to  substantiate  conclusions 
about  the  operation  of  the  INS.  Individual  flight  data  should 
be  in  a data  appendix  to  the  report  unless  specific  data  are 
required  to  substantiate  a conclusion  made  in  the  body  of  the 
report . 
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DR  NUMBER 


DATE 


VEHICLE  TYPE 

VEHICLE  SERIAL  KO(S). 

TEST  LOCATION 

MAJOR  SYSTEM/WU  C 

SU  B SYSTEM/WU  C 

DEFICIENCY 


DEFICIENCY  CIRCUMSTANCES /DESCRIPTION/CAUSES  (Continue  on  separate  page  it  naceaaary. ) 
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Figure  39  AFFTC  Deficiency  Report  Form 
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DATA  ANALYSIS 
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TERMINAL  ERROR  PLOT 

The  terminal  error  plot  shows  the  performance  of  an 
inertial  system  during  the  test  program  in  terms  of  shut- 
down error  for  all  flights.  A sample  plot  is  shown  in  Figure 
41.  Each  shutdown  error  obtained  during  the  test  program  is 
normalized  to  the  time  specified  in  the  INS  specification 
and  plotted.  The  CEP  is  calculated  and  plotted  as  a circle 
of  equal  probability  for  each  alignment  mode  evaluated.  The 
title  block  identifies  the  alignment  mode  and  the  time 
period  represented  on  the  plot.  The  specification  CEP  is 
shown  as  a circle  of  radius  equal  to  the  specified  value. 

The  data  for  the  terminal  error  plot  is  obtained  from 
the  flight  log  (Figure  40) . The  normalized  values  of  lati- 
tude error,  longitude  error  and  radial  error  are  calculated 
by  the  systems  test  engineer  and  stored  on  a computer  data 
card.  If  the  test  engineer  judges  that  the  data  from  a 
particular  flight  is  invalid,  it  should  be  documented  in  the 
inertial  platform  history  log  and  flight  log  and  not  punched 
onto  data  cards. 

When  a reasonable  amount  of  terminal  error  data  is  avail- 
able, the  data  is  run  through  a CEP  calculation.  After  the 
calculation  is  complete,  the  computer  examines  the  distribu- 
tion of  the  terminal  error  points  with  respect  to  the  value 
of  the  calculated  CEP.  If  the  terminal  error  from  any  of  the 
flights  exceed  a three-sigma  value  of  the  CEP,  that  flight 
terminal  error  is  temporarly  suppressed  from  the  calculation 
and  the  calculation  is  done  over.  The  test  engineer  should 
include  all  valid  terminal  error  data  that  is  available  in 
the  CEP  calculations.  New  flight  data  should  be  incorporated 
as  the  inertial  evaluation  portion  of  the  testing  progresses. 

POSITION  ERROR  PLOT 

If  inflight  accuracy  data  are  available,  position  error 
plots  showing  the  INS  accuracy  as  a function  of  elapsed  time 
should  be  produced.  Figure  42  shows  a sample  position  error 
plot  produced  from  actual  test  data.  The  plot  displays  lati- 
tude error,  longitude  error,  and  radial  error  as  a function 
of  elapsed  time.  On  the  plot,  north  latitude  error  and  east 
longitude  error  are  defined  as  positive.  Radial  error  is 
always  positive. 

The  AFFTC  computer  program  that  produces  the  position  error 
plot  is  called  the  Navigation  Analysis  Program  (NAVAN)#  In- 


mt 


puts  to  the  computer  for  a particular  flight  are  a radar  track- 
ing tape  of  that  flight  from  the  Space  Positioning  Branch,  an 
engineering  units  tape  produced  from  the  aircraft  instrumenta- 
tion raw  data  tape  for  the  same  flight,  and  card  entered  data 
for  information  that  can't  be  obtained  from  the  other  two 
sources.  The  program  merges  the  data,  edits  the  data  accord- 
ing to  given  editing  specifications,  writes  a data  file  of  the 
selected  and  edited  parameters,  performs  computations,  and 
produces  plots  of  the  computations. 

The  first  step  in  producing  position  error  data  is  to  make 
a time  correlation  table  for  the  flight  that  is  to  be  plotted. 
This  is  done  by  comparing  the  time  of  the  tone  cutoff  on  the 
Visacorder  data  recorded  at  the  Space  Positioning  facility 
during  flight  with  the  time  of  the  end  of  the  event  pulse 
from  the  onboard  data  for  each  test  point.  If  the  time 
difference  between  the  two  data  sources  shifts  significantly 
(greater  than  ten  percent  of  the  specification  drift  when  con- 
verted to  distance)  between  the  first  and  last  tone  breaks, 
the  data  will  be  unusable  unless  each  tone-break  is  correlated 
individually.  Both  data  sources  are  normalized  to  make  time 
equal  zero  be  the  time  that  NAV  was  selected  on  the  INS  control 
panel . 

Next,  the  data  is  run  through  a smoothing  routine  and 
interpolated  to  make  the  data  samples  from  the  two  sources 
occur  at  the  same  time.  The  smoothed  and  interpolated  data 
is  then  compared  to  obtain  the  latitude  error  and  longitude 
error.  Radial  error  is  calculated  from  the  latitude  and 
longitude  errors.  A listing  and  the  plotting  tape  required 
to  produce  the  position  error  plot  shown  in  Figure  42  is  then 
created . 

CEP  PLOT  FROM  POSITION  ERROR 

A single  CEP  plot  that  summarizes  the  results  of  all  of 
the  position  error  plots  £or  each  alignment  mode  is  shown  in 
Figure  43.  The  data  used  tro  generate  the  CEP  curve  are 
obtained  from  all  of  the  radial  error  information  from  all 
of  the  individual  accuracy  flights.  The  85  percent  confi- 
dence limits  shown  in  Figure  43  make  the  prediction  that 
85  percent  of  all  future  navigation  errors  will  fall  be- 
tween the  two  curves.  The  specification  line  is  the  accuracy 
that  the  contractor  was  required  to  meet.  This  number  is 
usually  a straight-line  specification  for  a finite  time-limit. 
Three  nautical  miles  per  hour  for  the  first  three  hours  of  ^ 
operation  is  the  specification  that  is  illustrated. 
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APPENDIX  A 


NAVAN  Computer  Program 


INTRODUCTION/PROGRAM  DESCRIPTION 


The  Navigation  Analysis  Program  (NAVAN)  is  a software  package 
developed  to  provide  the  test  engineer  with  a method  for  determining 
the  performance  of  inertial  navigation  systems.  NAVAN  was  designed  to 
analyze  INS  data  obtained  from  flight  test.  NAVAN  can  accept  data 
from  Radar  and  ADAS  tapes  and/or  punched  cards.  After  at  least  three 
flights  are  merged,  the  mean,  median,  R50,  R90,  CEP  and  confidence  limits 
are  calculated  and  data  is  plotted. 

Data  editing  of  the  Radar  data  is  done  by  checking  radar  velocity 
against  the  maximum  value  selected  by  the  engineer  (default  m 900  ft/sec) . 
If  the  velocity  exceeds  this  value  the  data  associated  with  this  time 
is  replaced  with  the  data  at  the  last  inbound  time.  Replaced  data  is 
noted  on  the  listing. 


To  match  Radar  and  ADAS  data  at  the  exact  same  times,  a linear 
interpolation  is  made  of  the  ADAS  data. 

Data  from  the  different  sources  will  not  occur  at  the  same  sample 
rate.  Because  a uniform  rate  is  required  for  statistical  analysis,  a 
least  squares  fit  is  made  of  the  data  and  values  from  all  sources  are 
established  at  the  same  sample  rate. 

The  initial  flight's  data  are  placed  on  the  New  Flight  File,  if 
not  the  initial  flight,  the  data  are  added  to  previous  data  and  are 
placed  on  the  New  History  Flight  File  (Old  History  Flight  File  on  next 
flight) . While  a time  history  plot  and  listing  of  the  data  are  avail- 
able, it  is  recommended  that  at  least  three  flights  be  obtained  before 
a statistical  analysis  is  made. 

THEORY 

The  methods  to  calculate  the  50th  and  90th  percentile  and  confidence 
limits  on  the  50th  percentile  used  in  the  analysis  program  are  presented 
in  this  section.  Two  methods  are  incorporated  in  the  program  for 
calculating  the  50th  percentile  value. 

The  following  definitions  are  used  throughout  this  section: 

x ■ latitude  error  (nm) 

y ■ longitude  error  (nm) 

r ■ radial  error  (✓ x 1 + ya)  (nm) 

m ■ number  of  testa 

i subscript  - test,  i ■ 1 thru  m 

q ■ a variable 

Eq^  ■ sum  of  a qj,  for  i ■ 1 thru  m 


Uq  = mean  or  expected  value  of  q 

a2g  = variance  of  q 

q = sample  mean  of  q = Zq-j/m 

S2g  = sample  variance  of  q = {Kq.^  - q)2}/m 

The  first  method  used  to  calculate  percentiles  of  radial  error  (R50 
and  R90)  is  based  on  the  Air  Standard  53/11. Bf  15  August  1968,  The 
Specification  and  Evaluation  of  the  Accuracy  of  Inertial  Navigation 
System.  The  computations  for  R50  and  R90  are  made  at  each  time  point  in 
the  following  manner: 

Suppose,  at  some  point  in  time,  there  are  rti  radial  errors  (r^,  1,  m) 
from  the  m corresponding  flights  in  the  sample. 

1.  Calculate  the  geometric  mean^(GM)  of  the  radial  errors: 


GM  = /Jlrj 


where:  H = the  product 


2.  Calculate  the  root  mean  square  i RMS)  of  the  radial  errors: 

RMS  = /Zr2i/m 

3.  Calculate  the  ratio  = GM/RMS 

4.  Define  RATIO  = GM/RMS  and  calculate  the  R50  and  R90  values  from: 


R50  a RMS 
R90  a RMS 


( . 7 RATIO  + .3) 

(1  + /I  - RATIO)  J 


R50  a RMS  ( . 7 RATIO  + . 4 /RATIO)  *\ 

R90  a RMS  {RATIO  + 1.6  (1  - RATIO2)}  J 


for  RATIO  < .6 


for  RATIO  > .6 


The  second  method  used  to  calculate  the  percentile  of  radial  error 
(CEP  only  in  this  program)  is  based  on  a paper  by  L.  L.  Rosen  and  D.  L. 
Harmer  titled  "Inertial  System  Performance  Evaluation"  which  was 
presented  at  the  Third  Inertial  Guidance  Test  Symposium  at  Holloman  AFB , 
New  Mexico,  in  1966.  A procedure  for  calculating  confidence  limits  on 
this  CEP  is  also  given  in  this  paper. 


from: 


At  each  time  point  the  percentiles  of  radial  error  are  calculated 


Rp  = °y  /a  (z'o2  + uz) 


Rp  = the  pth  percentile  of  radial  error  and  Zp  = the  pth  percentile 
point  of  a zero  mean  normal  distribution. 

Oy,  a,  oz,  yz  are  calculated  from  the  following  set  of  formulas: 

ox  « Sx/m/(m  - 1)  ; oy  a Sy/m/(m  - 1) 

lix  “ X ; yy  « y 

^Defined  in  glossary 
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K = ox/o y ; d = /mx2  + iiy1 

n = K2(2  - K2)  + 1 + (2/oy=)  (d*  jixlK2  - yyl) 

X = K2(K2  - 1)  + (l/oy2) (2yx2K2  + 2yy2  - d2) 
a = n + X ; b = X/a 

yz  = 1 - (2/9) (1  + b)/a  - (40/81) (b2/al) 
oz  - /(2/9) (1  + b)/a  + (16/27) (b2/a2) 

The  50th  percentile  (CEP)  is  calculated  from: 

Rp  = CEP  « Oy/aji*z 
since: 

Zp  = o for  p * 50 

At  each  time  point  the  100  (1  - o) % percent  confidence  limits  on  CEP 
(L1<CEP<I<2)  are  approximated  by  substituting  the  upper  and  lower 
confidence  limits  of  the  means  and  sigmas  calculated  below  into  the 
formulas  above : 

The  confidence  limits  on  the  means  are  calculated  using  a "t  test”:1 

|x|  - { (t<j/2)  (Ox//m)  } < ux  < | x | + {(to/2)  (ox//m)  ) 

|y|  - { (t0/2) (Oy//m) ) < My  < | x | + { (t0/2) (Oy//m) > 

The  confidence  limits  on  the  sigmas  are  calculated  using  a "chi 
test".1 

{Sx/m/(m  - i')'}/(Xi  - a/2)  < ox  < {Sx/m/(m  - l) }/xa/2 
{Sy/n/(m  - 1)  >/ ( Xi  - a/2)  < ay  < {Sy/m/(m  - 1)  }/xa/2 
where : 

ta/2  is  the  upper  100  (a/2)  percent  point  of  the  "t  distribution” 
with  m - 1 degrees  of  freedom.  Xa/2  and  “ a/ 2 are  the  lower  and  the 
upper  100  (a/2)  percent  points  of  the  "chiidistribution"  with  m - 1 
degrees  of  freedom. 

The  lower  limit  of  the  CEP  (Li)  is  calculated  by  using  the  lower 
limits  of  mean  and  sigma  to  calculate  = Oy/auz* . Should  the  lower 
limits  of  the  means  be  negative/  a zero  is  substituted  for  this  lower 
limit  in  the  computations  of  . 

The  upper  limit  of  the  CEP  (L2)  is  calculated  by  using  the  upper 
limits  of  mean  and  sigma  to  calculate  1*2  “ Oy/ayz* . 

PREPARATION  FOR  USE 

The  program  deck  of  NAVAN  will  be  permanently  stored  on  magnetic 
tape  number  04295  at  the  AFFTC  tape  library.  Prior  to  a first  use,  the 
compiled  file  must  be  copied  to  disc  and  stored  as  a permanent  file  as 
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shown  in  figure'  Al.  Once  the  program  is  stored  as  a permanent  file,  it 
may  be  attached  and  executed. 

SCOPD  OP  PROGRAM 

This  program  uses  control  information  from  data  cards  and  extracts 
selected  test  data  information  from  one  or  certain  combinations  of  the 
following  inputs:  (1)  An  Automatic  Data  Acquisition  System  (ADAS)  tape 

(produced  at  the  AFFTC) , (2)  A radar  tape  (produced  at  the  AFFTC) , (3) 

Card  data  containing  navigation  system  time,  latitude,  and  longitude,  and 
(4)  Card  system  data  containing  time,  latitude,  and  longitude  position 
error.  The  program  then  merges  the  data,  edits  wild  points,  performs 
various  calculations,  produces  plots  on  the  Calcomp  plotter,  and  produces 
printed  output  of  the  computations. 

Generalized  flowcharts  of  the  three  data  sources  (ADAS  - radar,  radar 
card,  and  system  card)  are  shown  in  figures  A2,  A4,  and  A6  respectively. 
Detailed  flowcharts  are  shown  in  figure  A29.  The  subroutines,  ADADATA, 
RADRD,  and  CRDDATA , select  the  times  and  parameters  from  the  data  sources. 
Multiple  events  on  the  ADAS  tape  are  read  with  one  start  and  stop  time. 

A wild  point  check  is  made  on  the  radar  data  only,  based  on  a maximum 
aircraft  horizontal  velocity  of  900  feet  per  second.  (Values  replaced 
are  noted  on  output  listing.)  This  value  may  be  changed  by  the  user  if 
required.  Any  radar  data  input  at  the  time  that  the  velocity  exceeds 
900  feet  per  second  is  replaced  with  the  last  inbound  value.  The  radar 
data  is  then  merged  with  the  ADAS  or  system  card  data,  based  on  the  time 
of  the  radar  data  point.  A linear  interpolation  is  used  to  obtain  the 
data  value  at  the  time  of  the  radar  data  point.  The  radar  latitude  and 
longitude  is  subtracted  from  the  system  latitude  and  longitude.  The 
resulting  error  is  converted  to  nautical  miles  and  the  radial  error  is 
computed  for  each  data  point.  A least  squares  curve  fit  is  then  made 
of  the  data  in  order  to  obtain  a constant  sample  rate.  The  sample  rate 
output  is  a user  defined  option. 

Note:  For  valid  statistical  data,  the  sample  rate  should  be  the  same 

for  all  flights,  regardless  of  the  data  source  (ADAS  - radar,  radar  card, 
or  system  card) . This  data  is  now  placed  on  the  New  Flight  File  (tape  16) . 
Tape  16  may  be  optionally  printed  or  plotted.  Tape  16  is  normally  the 
input  to  NAVMR.  This  is  combined  with  the  file  containing  data  from  the 
other  flights  to  be  analyzed  (tape  17) . This  merged  data  is  placed  on 
the  New  History  Flight  File  (tape  11),  and  may  also  be  optionally 
printed.  Tape  11  is  the  input  to  NAVAN  which  performs  the  statistical 
calculations  of  the  data  to  be  analyzed.  An  important  user  defined 
variable  in  NAVAN  is  the  time  span  over  which  each  calculation  occurs. 

This  time  span  should  be  equal  to  the  time  between  samples.  This  is  done 
to  obtain  only  one  data  value  per  flight  number  per  time  span,  in  order 
to  maintain  statistical  validity  for  the  calculated  confidence  intervals. 
NAVAN  generates  both  Calcomp  plots  of  the  data  and  a computer  listing. 


USER'S  GUIDE 


Listed  below  cure  the  necessary  tapes  for  options  being  used: 

LIST  OF  FILES 


FILE  NO.  DESCRIPTION 


4 Search  File  (Internal  Operating  File) 

5 Input  (Ccrd  Input) 

6 Output  (Printer  Output) 

7 Input  (ADAS  Tape) 

8 Input  (Radar  Tape) 

9 I/O  (ADAS  Selected  Parameters) 

10  Output  (Optional,  from  NAVAN  for  Debugging) 

11  New  History  Flight  File  (NHF,  New  Data  Base) 

12  I/O  (Statistical  Data  for  Plot,  Internal  Working 
File) 

13  Calcomp  Plot  Tape 

16  I/O  (New  Flight  File,  NFF) 

17  Input  (Old  History  File,  OHF,  Old  Data  Base) 

22  I/O  (Radar  Selected  Parameters) 

23  Internal  Working  File 


The  output  data  from  NAVAN  consists  of  an  output  file  computer 
listing,  as  in  figure  A8,  and  Calcomp  plots,  as  in  figures  A9,  A10,  and 
All.  The  New  History  Flight  File  (file  11)  contains  all  of  the  previously 
processed  flights  sequentially  arranged  in  order  of  increasing  time.  The 
output  listing  contains  a copy  of  all  of  the  input  cards  and  a copy  of 
all  of  the  processed  data.  This  can  be  used  to  check  the  validity  of  the 
data  prior  to  plotting. 


INPUT  CARDS  FOR  NAVAN 


The  cards  for  the  three  different  input  options  are  described 
separately. 


ADAS  - RADAR  COMBINATION  DATA 


NAVAN  data  card  #1,  figure  A12,  is  used  to  pick  the  input  option, 
INOP (2) «*1,  for  ADAS  and  Radar  tapes,  set  the  slope  of  the  CEP  line  to  be 
drawn  on  the  Calcomp  plots  and  select  other  program  options.  The 
NAMELIST  format  is  used  to  read  in  data.  A card  must  start  in  column  2 
with  $DATA  and  be  terminated  with  a $ sign.  Three  additional  input 
options  may  be  specified  on  this  card.  Normally  900  ft/sec  is  used 
internally  for  wild  point  checking,  but  this  may  be  reset  to  another 
value  by  placing  N900-XXXX.XX  (the  new  value)  on  the  card.  (Care  should 
be  taken  that  the  actual  aircraft  velocity  does  not  exceed  900  ft/sec.) 

NAVAN  data  card  <2,  figure  A13,  contains  the  number  of  parameters  to 
be  selected  from  the  ADAS  tape  and  the  parameter  ID  codes.  Note:  One 

data  word  may  be  made-up  of  two  words  containing  the  most  and  least 
significant  data  bits. 

NAVAN  data  card  #3,  figure  A14,  contains  data  to  correct  the  time  on 
the  ADAS  tape  to  the  time  on  the  radar  tape.  This  card  also  contains  the 
navigation  system  start  time. 


NAVAN  data  card  #4,  figure  A15,  contains  the  start  and  stop  search 
times  in  total  seconds. 

NAVAN  data  card  #5,  figure  A16,  is  in  NAMELIST  format  and  starts  with 
a $NAM1  in  column  2 and  is  terminated  with  a $ sign.  The  next  entry  is 
the  start  time,  in  seconds,  on  the  New  History  Flight  File,  tape  11,  where 
data  calculations  start.  Next,  the  stop  time,  in  seconds,  on  tape  11  is 
listed.  The  time  span  for  calculations  is  specified  and  should  be 
selected  so  as  to  have  one  data  point  per  time  span  per  flight  for  the 
best  statistical  validity  (standard  time  300  seconds) . 

NAVAN  data  cards  #6,  #6A,  #7,  #8,  and  #9,  figures  A17,  A18,  A19,  and 
A20,  are  only  required  for  Calcomp  plots.  Data  card  #6  has  the  X-axis  and 
Y-axis  scale  factors  for  each  plot,  the  initial  value  for  the  X-axis,  the 
number  of  plots,  and  a plot  heading  code.  The  number  of  plots  is  equal  to 
three  times  the  number  of  runs.  The  plot  heading  code  must  be  equal  to  1 
for  a heading  and  0 or  blank  for  no  heading.  Data  card  #6A  contains  the 
heading.  Data  card  #7  has  the  axes  lengths  for  Calcomp  plot  #1,  figure  A9, 
and  the  starting  value  for  the  Y-axis.  Data  card  #8  is  the  same  as  data 
card  #7,  but  contains  information  for  Calcomp  plot  #2,  figure  A10.  Data 
card  #9  is  the  same  as  data  cards  #7  and  #8,  but  contains  information  for 
Calcomp  plot  #3,  figure  All. 

If  there  is  to  be  more  than  one  data  run,  cards  with  the  new  values 
should  be  repeated.  The  end  of  the  data  request  should  be  indicated  by 
a blank  data  card,  figure  A21. 

SYSTEM  CARD  DATA 

This  is  the  case  where  system  error  data  will  be  entered  on  cards. 

NAVAN  data  card  #1,  figure  A22,  starts  with  a $DATA  in  column  2 and 
terminates  with  a $ sign.  The  NAMELIST  format  is  used  to  read  in  data. 

This  card  should  be  coded  with  INOP(3)=l  and  INOP(10)=1,  to  indicate  that 
the  data  cards  will  have  latitude  and  longitude  in  degrees.  Other  input 
options  are  the  same  as  for  ADAS  - Radar  Combination  Data,  figure  A12. 

NAVAN  data  card  #2,  figure  A23,  will  be  the  first  data  card.  One 
card  will  be  read  for  each  data  point,  and  the  final  data  value  will  be 
-100000.00,  figure  A24.  The  following  cards  after  the  final  data  value 
will  be  the  same  as  for  ADAS  - Radar  Combination  Data. 

RADAR  CARD  DAT* 

This  is  the  case  where  radar  tracking  data  are  available  and 
system  data  are  on  cards. 

NAVAN  data  card  #1,  figure  A25,  starts  with  a $DATA  in  column  2 and 
terminates  with  a $ sign.  The  NAMELIST  format  is  used  to  read  in  data. 

This  card  will  contain  the  codes  INOP(4)=l  and  INOP(10)=1,  with  the  other 
input  options  the  same  as  for  ADAS  - Radar  Combination  Data,  Figure  A12. 

NAVAN  data  card  #2,  figure  A26,  contains  the  start  and  stop  times 
for  the  radar  tape,  and  the  navigation  system  start  time. 

NAVAN  data  card  #3,  figure  A27,  is  the  first  data  card  and  should  be 
repeated  for  each  data  point,  and  the  final  data  value  will  be  -100000.00, 
figure  A28.  The  following  cards  after  the  final  data  value  will  be  the 
same  as  for  ADAS  - Radar  Combination  Data. 
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SAMPLE  DECK  SETUP  TO  RUN  ADAS-RADAR  COMBINATION  DATA. 

n 


77 


m 


1 FA K LAT 
i ED  IAN  (t 

fO 

<r<;^A  lat 
P«-S 

P 90 

Mt  AN  LUNG 
6"/RNS 
CEP 

SIGhA  LI 
P( 

90TH? 

■ *01*  CO 

NO.  DATA 

POINTS*  R 

, 76  25  00  F4  00 

. 601  5 11  E*  03 

-.2062  501*00 

.66  86  56 

0163121*00 

• TORI  736*01 

. R2 13  60  E*0C 

.72  72  86 

98  **C3^(* 

. o^^Aij  fflNP  TtVf - 

. 133388E+01 

.10R2RRE*C1 
»-«  (L-ARfc r 38608 

.176068 
3E  *001  LOWE 

rt  EA  N R 
P(2> 


SIGNA  K 
P(  3) 


P<A) 


• 100711E+01 
• 10R660E+01 


• 62 5303E+00 
• 1R0312E*01 


.1752176*01 


PERCENT  CONFIDENCE  LllUTS  HR  CEP  ARt 
PERCENT  Cn  NF  I D FN  CF  LINITS  FIR  CEP  ARt 
. 1086tlf*01  0.  0. 


,AlA610E*0u<  LJWtr.1 
. R65165E*0L'<L3Wt  PI 

0. 


AND 

AND 

AND 


.A53326E*01<UPPEP) 
.35eR31E*01tUPPER  1 
,270121E*01(UPPER> 
0.  0. 


0. 


PERCENT  CPNFIDFNCr  LINITS  FOR  P70TH  ARE  . 70C 32&E +OOANO  .B03132E*ul 


M J» 

l T* 

I G* 

p--*_ 


.ACOOOOf  *01 
.30000  OF  *01 
■.762  500F  *00 
-.7062501*00 


l A» 

I.C* 


l A* 

LX* 


• 1 OR  175E  *0 1 
.108  5A  51*01 

• 1 35  IT  OE  *0 1 
.T  17  30  IE  *07 
.001  81  IE  *00 


.‘25  30  3e  *00 
• 1 OR  ?A  AF  *0 1 
.°10  31  ?E  *0  0 


Figure  A8  Program  navan  Computer  Output  Listing 
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Input  Options* 


INOP (2) =1 
INOP  ( 3)  **1 
INOP (4) *1 
INOP (5) -1 


INOP (6) =1 


INOP (7) =1 


ADAS  - Radar  Combination  Data 

System  Card  Data 

Radar  Card  Data 

Prints  ADAS  and  Radar  Data 

before  merging 

Prints  data  written  onto 

tape  16  (NFF) 

Produces  Calcomp  plot  of  data 
on  tape  16  (Calcomp  plot  #1) 


Slope  of  CEP  Specification  line  for 
Calcomp  plot 


SLOPE *XX . XXXX  (nm/unit  time  of  plot  absissa) 


N900=XXXX.XXXXX  (feet/second) 


A new  value  entered  here  changes  the  value 
of  velocity  used  for  wild  point  editing  of 
radar  data. 


Termination  of  card  by  a $ sign. 


* This  card  is  of  the  NAMELIST  format. 


Figure  A12  Program  NAVAN  Data  Card  #1 

•2 


I " 


nifiaannontrnnnnnnann-innnannanaTnannnnnnnaiinnnnaannTnnBnnaonnrrncanaa 


IE 


nPDonnPDQl 


(F10. 3) 


Stop  time  (total  seconds) 


(F10.3) 


$NAM1 


Beginning  time  (seconds)  on  New  History 
Flight  File  (tape  11) 


BEG-XXXX.XXXX 


Stop  time  (seconds)  on  New  History  Flight 
File  (tape  11) 


STP=XXXX.XXXX 


Time  interval  between  data  points  (seconds) 


DTIME=XXXXX . XXX 


Termination  of  card  by  a $ sign. 


NOTE:  This  card  is  of  the  NAMELIST  format. 


Figure  A16 


Program  NAVAN  Data  Card  #5 
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Figure  A17 


Program  NAVAN  Data  Card  #6  and  #6A 


X-axis  length  (inches)  for  plot  no.  1 


(latitude  error,  longitude  error,  and 
radial  error)  (F10 . 3) 


Y-axis  length  (inches)  for  plot  no.  1 


Starting  value  of  Y-axis  for  plot  no.  1 


COMMENT:  Starting  value  for  X-axis  is  the 

same  for  all  3 plots. 


Figure  A18  Program  NAVAN  Data  Card  #7 
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Bj  Input  Options* 
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If  blank: 
INOP(10)=1 

INOP (4) =1 


data  is  in  nautical  miles, 
latitude  and  longitude  data 
are  in  degrees. 

Radar  Card  Data 


Slope  of  CEP  Specification  line  for 
Calcomp  plot 

SLOPEsXX . XXXX  (nm/unit  time  of  plot  absissa) 


Termination  of  card  by  a $ sign. 
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NOTE:  This  card  is  of  the  NAMELIST  format. | 

* Input  Options  are  the  same  as  for  Program 
NAVAN  Data  Card  #1. 
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Figure  A25 


Prograun  NAVAN  Data  Card  11 


Elapsed  time  (seconds) 


NOTE:  Repeat  this  data  card  for  each  data 

point. 
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Program  NAVAN  Flowchart 
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Figure  A29 


Program  NAVAN  Flowchart  (continued) 

IN 


SUBROUTINE  XPLOTS 


Figure  A29  Program  NAVAN  Flowchart  (continued) 
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Prograa  NAVAN  Flowchart  (continued) 
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CEPLOT  Computer  Program 


III 


INTRODUCTION 


The  Circular  Error  of  Probability  Program  (CEPLOT)  calculates  the 
CEP  and  plots  the  data.  The  input  is  navigation  system  end  point  position 
error  data  normalized  to  one  hour,  based  on  a linear  error  rate.  Data  are 
input  in  card  form. 

OVERVIEW 

CEPLOT  consists  of  the  main  program  and  subprograms . It  reads  the 
data,  computes  the  error,  normalizes  the  data,  computes  the  CEP,  and 
produces  a tabular  listing  and  Calcomp  plot  tape  for  plotting. 

A flowchart,  is  shown  in  figure  Bl- 

The  method  used  to  calculate  the  CEP  is  the  same  as  is  used  in  NAVAN 
(Appendix  A),  which  is  based  on  the  Air  Standard  53/11. B,  15  August  1968, 
The  Specification  and  Evaluation  of  the  Accuracy  of  Inertial  Navigation 
System. 

PREPARATION  FOR  USE 

The  program  deck  of  CEPLOT  will  be  permanently  stored  at  the  AFFTC 
Systems  Engineering  Branch. 

USER'S  GUIDE 

The  deck  developed  for  running  CEPLOT  is  composed  of  three  parts. 

The  first  part  is  the  job  control  cards  which  are  used  to  attach  the 
Calcomp  plot  tapes  and  execute  the  program.  Examples  of  these  control 
cards  are  shown  in  figure  B2. 

The  second  portion  of  the  program  is  the  main  program  and 
subroutines . 

The  third  portion  of  the  program  consists  of  two  parts:  (1)  Tne 

parameter  cards  which  control  the  program  and  plotting,  and  (2)  The  data 
cards.  These  data  cards  are  shown  in  figure  B3. 

The  output  data  from  CEPLOT  consists  of  a computer  output  listing, 
figure  B4,  and  a Calcomp  plot,  figure  B5. 

INPUT  CARDS  FOR  CEPLOT 

CEPLOT  data  card  #1,  figure  B6,  contains  the  number  of  data  points  to 
be  calculated  and  plotted,  and  the  heading  to  be  placed  on  the  Calcomp 
plot. 


CEPLOT  data  card  #2,  figure  B7,  contains  the  type  of  INS  alignment, 
either  gyrocompass  (GC) , best  azimuth  heading  (BATH) , or  standard  (STD) , 
and  the  format  of  the  input  data.  The  data  format  may  be  "RAW"  which 
means  that  the  delta  latitude,  delta  longitude,  flight  number,  and  time 
will  be  input.  Radial  error  will  be  calculated  and  the  data  normalized. 

In  the  other  format,  the  data  are  input  as  delta  latitude,  delta  longitude, 
and  radial  error  (PMT  Left  Blank) . 

CEPLOT  data  card  #3,  figure  B8,  the  first  data  card,  if  not  "RAW", 
will  contain  delta  latitude,  delta  longitude,  and  radial  error  which  has 


Ml 


been  normalized.  One  card  is  used  for  each  data  point.  If  "RAW"  data, 
the  card,  figure  B9,  will  have  flight  number,  delta  latitude,  delta 
longitude,  and  time  in  hours  and  minutes.  This  data  will  not  have  been 
normalized.  One  card  is  required  for  each  data  point. 


FLOWCHART  OF  PROGRAM  CEPLOT 


FLOWCHART  OF  PROGRAM  CEPLOT 
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FLOWCHART  OF  SUBROUTINE  BPLOT 


Figure  B1  Program  CEPLOT  Flowchart  (concluded) 
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Figure  B5  Program  CEPLCtT  Caloonp  OufcputvFlot  ,„r.  i? 
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Figure  B7b'i/  Program  CEPLOT  Data  Card  #2 
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Flight  number 


Delta  latitude  (degrees) 


(F9.3) 


Delta  longitude  (degrees) 


Time  (hours) 


Time  (minutes) 
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Figure  B9  Program  CEPLOT  Data  Card  #3  and 
Following  Cards 
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APPENDIX  C 
Check  Cases 


The  following  check  case  for  program  NAVAN  uses  data  from  an 
onboard  navigation  system  and  space  positioning  (RADAR)  tracking  data. 
The  ADAS  and  RADAR  data  are  first  merged  and  placed  on  the  New  History 
Flight  File,  then  another  set  of  data  is  added  to  this  file  and  the 
merged  data  is  analyzed  and  the  results  are  plotted.  All  intermediate 
data  was  printed  for  help  in  diagnosis  of  program  errors. 
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The  check  case  for  program  CEPLOT  is  a set  of  random  end  point 
errors  (generated  using  a random  number  routine)  and  illustrates  the 
program  operation  including  plotting  of  the  data. 
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RAW 

01 -.*9  5 
£12.55  5 
01 -.  95  1 
011.399 
£11.02 
01 1.32  A 
911.A97  - 

feft 

01.5  75 

% -1806 
01 -.32  A 

• A 39 

• 1 A3 
_ -.  71 7 

0i-l  .9  55 
01 -.53  6 
OH. 65u 
*i  l.Oi  * - 

“1-1.3  37 

<n.i9c 

Si  IS 1 

bi:iu 

01-1.019 
END  MA  5S 
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PROGRAM  CEPLOT  CHECK  CASE  INPUT 


GLOSSARY 


anisoelasticity 


azimuth 
azimuth  angle 


chi  distribution 
test 

circle  of  equal 
probability  (CEP) 


cross  coupling 


geometric  mean 


heading 


elastic  deformation  of  a gyro's  gimbals  that  result 
in  gyro  drift  proportional  to  the  g2  value  of 
vibration. 

horizontal  direction  or  bearing. 

azimuth  measured  from  0°  at  the  north  or  south 
reference  direction  clockwise  or  counterclockwise 
through  90°  or  180°. 

a statistical  significance  test  based  on  frequency 
of  occurance. 

a measure  of  the  accuracy  with  which  an  aircraft 
can  be  guided;  the  radius  of  the  circle  at  a 
specific  distance  in  which  50  percent  of  navigation 
errors  fall;  also  called  circular  error  probable, 
and  circle  of  probable  error. 

The  position  of  the  accelerometer  pendulum  is  off 
null  by  a small  angle  in  order  to  create  the  error 
signal  necessary  to  rebalance  the  accelerometer. 
This  angle  causes  cross-coupling  of  the  component 
of  acceleration  which  is  normal  to  the 
accelerometer . 

a measure  of  central  position.  The  geometric  mean 
of  n quantities  equals  the  nth  root  of  the  product 
of  the  quantities. 

horizontal  direction  in  which  an  aircraft  is 
pointed,  expressed  as  angular  distance  from  a 
reference  direction. 


inertial  coordinate 
system 


latitude 


a system  in  which  the  (vector)  momentum  of  a 
particle  is  conserved  in  the  absence  of  external 
forces.  Thus,  only  in  an  inertial  system  can 
Newton's  Laws  of  Motion  be  appropriately  applied. 

Terrestrial  latitude  is  angular  distance  from  the 
Equator,  measured  northward  or  southward  through 
90°. 


longitude  Terrestrial  longitude  is  the  arc  of  a parallel,  or 

the  angle  at  the  pole,  between  the  prime  meridian 
and  the  meridian  of  a point  on  the  Earth,  measured 
eastward  or  westward  from  the  prime  meridian 
through  180°. 


in 


¥ 


IP 


X 


magnetic  pole  either  of  the  two  places  on  the  surface  of  the 

Earth  where  the  magnetic  dip  is  90° , that  in  the 
Northern  Hemisphere  (at,  approximately,  latitude, 
73  8 N,  longitude,  101  W in  1955)  being  designated 
North  Magnetic  Pole,  and  that  in  the  Southern 
Hemisphere  (at.  approximately,  latitude,  68°S, 
longitude,  144  E in  1955)  being  designated  South 
Magnetic  Pole. 

proportional  bias  the  bias  required  to  prevent  gyro  drift  that 

results  from  acceleration  perpendicular  to  the 
output  axis  at  the  gyro. 

root  mean  square  in  statistics,  the  square  root  of  the  arithmetic 

error  (RMS)  mean  of  the  squares  of  the  deviations  of  the 

various  items  from  the  arithmetic  mean  of  the 
whole,  also  termed  standard  deviation. 

Student's  t a statistical  method  to  test  the  hypothesis  that 

distribution  test  the  mean  of  the  sample  is  consistant  with  being 

equal  to  the  mean,  y,  of  the  assumed  population. 

vibropendulous  the  accelerometer  error  produced  when  the 

pendulous  mass  in  an  accelerometer  is  vibrated. 
This  error  is  proportional  to  the  vibropendulous 
coefficient  times  this  vibration  in  feet  per 
second  per  second. 
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LIFT  OF  ABBREVIATIONS  AND  SYMBOLS 


Item 

Definition 

Unit 

a 

( *▼  T^Ot>  Vj  XTi) 

acceleration 

ft  per  sec2 

ac 

centripetal  acceleration 

ft  per  sec2 

d 

distance 

NM  or  ft 

E 

East 

- - - 

g 

acceleration  of  gravity  at  mean  sea  level 

32.3  ft  per  sec 

GM 

geometric  mean 

dimensionless 

i 

(subscript)  test,  i = 1,  m 

dimensionless 

L 

confidence  limits 

dimensionless 

m 

number  of  tests 

dimensionless 

II 

North 

- - - 

q 

a variable 

dimensionless 

q : 

sample  mean  of  q 

dimensionless 

r 

radius  of  the  Earth 

2.09  x 107  ft 

r 

radial  error  = V x2  + y2 

NM 

RATIO 

GM/RMS 

dimensionless 

R 

P 

pth  percentile  of  radial  error 

dimensionless 

S 

South 

- - - 

t 

time 

•thy  t 

sec 

t 

t in  Student’s  t test 

- - - 

T.O. 

Technical  Order 

- - - 

V 

velocity 

ft  per  sec 

W 

best 

- - - 

X 

latitude  error 

NM 

y 

lonoitude  error 

NM 

z 

vertical 

- - - 

z ' 
p 

pth  percentile  point  of  a zero  mean 
normal  distribution 

dimensionless 

—r 


Item 

Definition 

Unit 

o 

(l-o)%  confidence  intervals 

dimensionless 

*a 

azimuth  drift  rate 

deg  per  hr 

Ea 

acceleration  error 

ft  per  sec* 

ed 

distance  error 

NM  or  ft 

level  gyro  drift  rate 

deg  per  hr 

ev 

velocity  error  inside 

the  Schuler 

loop 

ft  per  sec 

evo 

velocity  error  outside 

the  Schuler 

loop 

ft  per  sec 

eo 

azimuth  angle  error 

deg 

ee 

pitch  angle  error 

deg 

% 

roll  angle  error 

deg 

e 

pitch  angle 

deg 

X* 

<j/r 

l/secl 

"q 

mean  or  expected  value 

of  q 

dimensionless 

s 

sample  sigma 

dimensionless 

V 

sample  variance  of  g 

dimensionless 

o 

azimuth  angle 

deg 

variance  of  q 

dimensionless 

sum  of  q^  for  i * 1,  m 

dimensionless 

♦ 

latitude 

deg 

<• 

roll  angle 

deg 

X 

chi 

dimensionless 

<i > 

gyro  torquing  rate 

deg  per  hr 

n 

Earth  rotation  rate 

15.04  deg  per  hr 

U) 

radians  per  unit  time 

rad  per  hr 

1M 


